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NOMENCLATURE 

All symbols are defined whenever they appear first 
in the text. 

a - length of plate in the current state 

a 0 - initial length of plate 

a^,a 2» a 5 - principal stretches 

b - width of plate in the current state 

b Q - initial width of plate 

c 11» c 12**“ moduli in eq.. (3.2) 

c 11t o 12 ..- moduli in eq.. (3.4) 

^ijkl “ instantaneous moduli in Chapter 2 

D - determinent 

E - Young's modulus 

e - compressive sti*ain 

ey - compressive strain at yield 

f( ) - probability density function (Appendix A) 

E|| - Mises critical stress function 

Erj, - Tresca critical stress function 

Ejj - analytical minimum stress function for Mises 

E* - analytical minimum stress function for Tresca 

E( ) - cumulative probability density function(Appendix A) 

0'1>&2» g '3 “ shear moduli 

h - rate of hardening 

h a p - hardening matrix 

J - stress invariant (= a^a 2 a^) 
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- probabilistic critical stress functions (Appendix A) 

- mean values of critical stress function (Appendix A) 

- moduli defined in eq. (3.20) 

- number of half-waves along lengthu at buckling 

- components of unit normal to yield surface 

- unit surface normal to the element 

- indeterminable stress parameter eq. (4.7 ) 

- probability of failure (Appendix A) 

- reliability 

- surface 

- nominal stress tensor 

- standard deviation of the variable (Appendix A) 

- thickness of the plate in the current state 

- initial thickness of the plate 

- surface traction 

- volume 

- coefficient of variation (Appendix A) 

- velocity components in x^ directions 

- Ogden's strain energy function 

- velocity components along x,y,z direction 

- coordinate system in suffix notation 

- term indicating probability of failure (Appendix A) 

- Ogden's parameter in strain energy f motion . 

- parameter in the stress-strain relation eq. (3.45) 

- shear rate 
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coupled hardening parameter in Chapter 3 
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shear modulus in Chapter 3 
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true (Cauchy) stress 
yield stress 

a parameter defined by eq_. (3.49) 
strain -rate tensor 
rotation tensor 
mb/a (Appendix A) 
material derivative 
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ABSTRACT 


The problem of buckling in thin rectangular plates 
under uniaxial compression is studied when (i) the material 
is elastic-plastic following the Tresca yield criterion 
with coupled hardening between the yield surface facets 
which meet at a vertex; (ii) the material is elastically 
incompressible and develops anisotropy during deformation. 
In the elastic-plastic case it is shown that the yield 
surface comer considerably lowers the bifurcation stress, 
either for independent hardening or for coupled hardening, 
compared to that obtained by using the Mises yield 
criterion; the latter is known to give buckling stresses 
much higher than the experimental results. In elastic 
case, the stress at bifurcation is found for a material 
defined by the Ogden's strain energy function, lastly, 
the inelastic buckling problem has been re-examined with 
a probablistic approach; the buckling stress has been 
calculated for a given probability of failure and for a 
desired variation in the material parameters and the 
plate dimensions. 



CHAPTER 1 


INTRODUCTION 

1 . 1 BRIEF LITERATURE REVIEW 

I"t is not the aim here to review a vast amount 

of literature available in the theory of stability of 
structures, mostly under ’small deformation’ approximation. 
Related to the present work we concentrate only on the 
subject matter based on the continuum approach to the 
stability investigation. 

A new trend in the subject of stability was 
induced due to Shanley’s (1947) demonstration of the 
fact that the problems of bifurcation (i.e. non-uniqueness 
of solution) and that of instability are, in general, 
two different concepts though there is a definite relation- 
ship between each other. Some ten years later, Hill 
presented a general three-dimensional continuum theory 
of uniqueness and stability, and thereby established 
exact connections between the uniqueness of deformation 
and the stability of equilibrium configuration for elastic, 
rigid-plastic and elastic -plastic solids. These principles 
are also available in a compact form (Hill, 1959, 61 ). 

In examining the problem of uniqueness of a 
body under finite strain, usually two methods have been 
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adopted. The first is the method of small deformations 
superimposed on finite deformations as available in, 
for example Green and Zema(l954) and Biot (1965). A 
system of equations is obtained for the equilibrium of 
the incremental deformation; these equations together 
with rate boundary conditions and rate material properties 
constitute the rate problem. The second method makes use 
of a variational technique in which a certain functional 
of velocity gradients (or of stress-rates) vanishes in 
the absence of a unique solution; non-vanishing of the 
functional implies uniqueness. A detailed account df 
both the approaches is contained in Hill f s work. 

Solutions of some specific, rate problems in 
elastic solids have been given by a number of investigators, 
most notably by Biot (1965). Plain strain/plane stress 
bifurcations in elastic-plastic plate were examined, 
amongst many, by Ramsey (1969), Dubey and Ariaratnam(l969) , 
Hill and Hutchinson (1975), Young (1976) and Needleman 
(1979). Wilkes (1955), and later on Vaughan ( 1 97 1 ) and 
Sierakowski et.al, (1975), obtained bifurcation conditions 
in solid elastic cylinders under end thrust, while insta- 
bilities in thin tubes were investigated by Dubey (1969) 
and Storakers (1975). Recent studies by Haughton and Ogden 
(1979a, b) and ITiyogi (1980) present a comprehensive treatment 
of instabilities in thin and thick-walled elastic cylinders 
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under axial compression. Bifurcations in elastic plastic 
plates under pure bending have recently been reported by 
Triantafyllidis (1980). 

Hill and Sewell (i960 a,b) have applied the 
variational technique in an investigation of the buckling 
of an inelastic column. The analysis has been extended 
to the buckling of elastic -pi as tic plates also (Sewell 
1963, 64 ). Some investigations in the stability of rigid- 
plastic cylindrical shells under pressure and/or axial 
load are available due to Stor5kers (1971), Kumar and 
Ariaratnam (1975,76) and Kumar and Shukla (1981). The 
finite element formulation of Hills variational principle 
has been successfully applied to study instabilities in 
cylindrical bars under tension and compression (Feedlemen 
1972, 73), bifurcation in plane tension test (Burke and 
Nix, 1979), bifurcation in cylindrical shell under 
internal pressure (Chu, 1979) and buckling of plates in 
the presence of cracks (Markstrom and Stor&kers, 1980). 

1.2 PRESENT WORK 

The thesis presents an application of Hill’s 
variational technique to examine stability in axially 
compressed thin elastic and elastic -plastic plates. 

There are large number of investigations available 
concerning buckling of plates in the inelastic range. 
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but most of them are -merely an extension of the classical 
buckling theory in elasticity. The most notable contri- 
butions are due to Sewell (1963,64) who applied the 
continuum theory of uniqueness and stability to determine 
the buckling stress of an elastic plastic plate of Mises 
material. It was also shown that buckling stress could 
be sensitive to changes in the directions of the normals 
at the Ire sea yield surface corner and that it could 
predict a value of the critical stress much lower than 
that for the Mises yield surface. Similar observations 
were also reported later on by Dubey and Ariaratnam(l969) . 
The present investigation reexamines Sewells problem vixen 
the current state of stress is at the comer of Tresea 
yield surface. The material is supposed to obey Hill's 
(1966) constitutive relations for metal crystals defor- 
ming by multislip. Then, using the flow-rule at a 
vertex (Sewell, 1972, 74), it has also been possible to 
include the effect of coupled hardening between the facts 
of the Tresea yield surface that meet at a vertex. It 
is shown, by using a stress-strain curve of the Ramberg- 
Osgood type (or even otherwise) that a considerable 
reduction in the critical stress is possible by using the 
Tresea yield criterion instead of the Mises, which is 
known to give higher buckling stress as compared to the 
experimental values. 
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A probablistic approach, in brief, is also 

cJL 

presentsgg to study the buckling of elastic-plastic 
plates. Curves are drawn to obtain the buckling stress 
for a known probability of failure and for the desired 
variation in the material parameters and the plate 
dimensions. 

Also, buckling of thin elastic plates has 
been reviewed by considering Ogden^s (1972) strain-energy 
function for a general incompressible material. The 
material is initially isotropic and develops anisotropy 
due to deformation. The outline of the thesis is as follows: 

Chapter two presents a brief account of Hill's 

* 

sufficient conditions for the uniqueness of deformation 
and the stability of equilibrium configuration. The 
stability criterion is applied to an axially compressed 
thin elastic -plastic plate in Chapter three and to an 
elastic plate in Chapter four. The conclusions are 
summarised in Chapter five. Appendix A gives a 
probabilistic approach to determine the buckling stress 
of an inelastic plate. Appendix B gives the standard 
normal probability tables. 



CHAPTER 2 


GENERAL THEORY OP STABILITY AND UNIQUENESS 

The process of continuing deformation of a 
finitely deformed solid is marked by large elastic or/ and 
plastic deformations terminating in non-uniqueness, Upto 
what extent the deformation of a structure is unique, under 
given loading conditions, can be determined by using some 
'uniqueness condition' and when the deformation becomes 
non-unique, i,e. when a bifurcation occurs, can be found 
by solving the rate boundary value problem. In this 
chapter, we shall discuss the formulation of the rate 
problem and derive sufficient conditions for uniqueness 
of deformation and stability of equilibrium configuration. 

2.1 RATE PROBLEM 

It is usually convenient to represent the equations 
of equilibrium for the incremental deformation from some 
current equilibrium state in the rate form. Similarly, the 
boundary conditions and the constitutive relations for the 
incremental deformation can also be expressed in the rate 
form. The rate equilibrium equations together with the 
rate boundary conditions and the rate material properties 
constitute the rate problem. 
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2. 1 . t E quilibrium Squat ions 


Set the body in an initial configuration B° be 
finitely deformed to some current configuration B # 
Suppose that in the current state, the traction T.= are 

J 

specified on some part of the surface S, i.e. 


1 5 = n i °13 (2-1) 

and the displacements are prescribed on the remainder. 
Here, in (2.1 ) <7^ is the true (Cauchy) stress and nj_ 

is the unit normal to the element. Then, if the body 
force is per unit volume in the current configuration, 
the equations of equilibrium are 


<y . . . + f . = o 

10,1 0 


( 2 . 2 ) 


where a comma denotes derivative with respect to ; 
or alternatively can be expressed as 

s ij,i + = 0 (2.2a) 

when the current state is taken as the reference state. 

In (2.2a) s^- is the nominal (lagrangian) stress. 

Bor the analysis of the problem of bifurcation 

of equilibrium, we need to examine the behaviour of 

the body when it is perturbed from its deformed state B 

to a neighbouring state B ' . During this incremental 

deformation, let the velocity v^ and the nominal traction 

rates , t i be prescribed on part Sy and S T respectively, of 
J 
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the surface S of the body* The equations of continuing 
equilibrium in terms of the lagrangian stress-rate are 
(Hill, 195 9) 


s - • = 0 in V 


(2.3) 


neglecting body forces. Here, 7 is the current volume 
of the body. The material derivative of the true stress, 
^i j -*- s related to by the following relation 


iD 


a 


ij v k,k 


- o. 


3 ^ 


v i,k 


(2.4) 


when the material is incompressible (i.e. v^^ = 0) and 
the current stress-distribution is homogeneous, the 
equilibrium equation can alternatively be expressed as 



i 


0 in V 


(2.3a) 


Equations (2.3) are the most concise form of the rate 
equilibrium relations. They have, however, been expressed 
in various other forms in the literature* Some of the 
well-known ones have been discussed by Bazant (1971). 

Bor example, in terms of Jaumann flux of the Oauchy 
stress, we have 


DO.. 

(""imT" ^ + ^ 'lk “jk + a ij v k,k” ^k e ik\i = 0 

(2.5) 

as is available in Biot (1965). Here (Vj_ .p V D*,i^ 

is the antisymmetric part of the velocity gradient tensor 
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while S^j = t (v.^ + is the symmetric part. On 

the other hand, Green and Zerna (1954-) and Wesolowski 
(1962) express their basic equations in terms of convected 
derivative of the Gauchy stress 


6 <? 


ID 


6t 


),i + ( 


a. .v- 


iO v k,k + 


ik v ;j,k\i “ 0 (2 * 6 ^ 


2.1,2 Boundary Conditions 


The rate boundary condition are normally written as 


T • 
X D 

= n i 

on S, 

v i 

= V i 

on S' 


The quantity v^ is usually a given function of position 
on the boundary, while T-j (the nominal traction-rate) 
may depend on position and local unknown velocities. 

Various formsfor T j , involving configuration dependent 
loading, have been derived by Hill (1962) and Dubey(l970) . 
For example, the dead loads are characterized by T-j = 0. 

2.1.3 Constitutive Relations 

We consider only those elastic-plastic solids having 
no natural time (i.e. no creep or real strain-rate effects). 
The inst ant aneous (incremental ) material behaviour ensuing 
from the given current state is expressed by relation of 
the form 

Stress-rate = f (strain-rate) (2.8) 



where f is a homogeneous function of degree one in the 
strain-rate components with the property f (0) = 0, but 
may be non-linear. At this stage we need not specify 
the forms of *f ’ which will be done at a later stage. 

2.2 UNIQUENESS CRITERION 


Eor incremental deformation from the current 
configuration B to a Neighbouring configuration B f , we 
have relations (2.3) and (2.7), A condition for unique- 
ness of solution is obtained as follows. Suppose that 
there could be more than one solution of a bound ary-value 
problem and let the difference between any two correspond- 
ing quantities be denoted by the prefix A . Then 


A S - - ■ = 0 in V 

J- J , X 

(2.9) 

A $ j = n i ^ ^ -j = 0 on S T 

(2.10) 

AT. =0 on s v 

(2.11) 


Hence 


/ 

S 



A v . dS 

J 



A s.- -jAv . . 
13 3,i 



( 2 . 12 ) 


by using Gauss's divergence theorem. A sufficient condition 
for uniqueness is, therefore, that 


f as. . 4y- . dV > 0 (2.13) 

Y 13 3,i 

for all pairs of continuous differentiable velocity fields 
taking the prescribed values on S and the corresponding 
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stress-rates s ± ^ obtained from the v ± through the 
constitutive relation. The reverse inequality would also 
give a sufficient condition, but is not considered here 
because such states are unstable (shown later on). It is 
important to note that when the constitutive relations 
are linear, the integral in (2.13) is a single-valued 
functional of v^ (regarded as a single field), but when 
they are non-linear, it is multivalued, 

2.3 STABILITY CRITERION 

To be definite suppose that a part S y of the 
surface of the body is rigidly constrained and that 
dead loads are maintained on the part Srj. The body, 
initially at rest, is imagined to be set in motion by 
some kind of perturbation . The initial state is said 
to be stable if the amplitude of the ensuing displace- 
ment is always vanishingly small whenever the perturba- 
tion itself is. If, on the other hand, the amplitude is 
finite for atleast one type of disturbance, however small 
this might be, the state is said to be unstable. Evidently, 
a sufficient condition for stability is that the internal 
energy stored should exceed the work done by the external 
loads during the resulting motion [of., e.g. Hill (1959)» 
Prager (1961)] , 

In a time interval 6t, due to a virtual displacement 
Vj_ 6t, the nominal stress changes from Sj_j to ( s ij +s ij6t), 
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still referred to the current state. Therefore, the 
increase in the internal energy will be 


6t v (s i3 + 1 6t) v j,i dV 

The virtual work done by the external forces 
during the time interval 6t, under the load on Srp and 
rigid constraints on Sy, i.e. Vj = 0 on , is 


6t f T.y.dS = 6t / n.s, i T i dS = 6t f s. -v • -dV 
S JJ s l 13 3 y ID 3»i 


using the divergence theorem and the current equilibrium 
condition = 0 . Hence, to the second order in St, 
the internal work minus the external work is 



(2.14) 


How, suppose that K 0 is the initial infinitesimal kinetic 
energy of the body. After the time St, in the new confi- 
guration, its energy is K = K Q - W, by the principle of 
conservation of energy. Since K can never be negative, 
an upper bound for St can be obtained using (2,14) if 
W remains positive. Therefore, the positive— definiteness 
of (2.14) is sufficient for stability, i.e. for stability 
we must have 


f 

Y 



v. , dV > 0 
3» 1 


(2.15) 


for all continuous differentiable velocity fields which 
vanish on S y (but are not identically zero) and are 
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compatible via. the constitutive equations with the 
stress rates 

2.3.1 Stability vs. Uniqueness 

Consider the relation between the conditions ( 2 . 1 3 ) 
and (2.15). We note that the Av i fields admitted in 
(2,13) are exactly those admitted in (2.15) since both 
vanish on S v and are otherwise arbitrary. However, (2.15) 
is a single-valued functional while ( 2 , 13 ) is multi-valued 
when constitutive equations are non-linear and single- 
valued only when they are linear . Consequently, the 
conditions are identical only when the rate -equations 
are linear. 

However, it may be mentioned here that the 
condition ( 2 . 13 ) is not necessary for uniqueness; the 
positive identification of bifurcation states involves 
actually demonstrating solutions to ( 2 . 9 ) - (2.1 l). For 
certain linear solids, a solution does exist in those 
states for which the functional in ( 2 . 13 ) is generally 
positive but vanishes for at least one non-zero field 

A v • • It is such bifurcations that would normally be 
J 

reached first in any actual process of deformations since 
they terminate a continuous range of states which not 
only satisfy ( 2 . 1 3 ) but are also stable. 

For the linear elastic— plastic solids, it is 
obvious now that the conditions for uniqueness and 
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stability are identical ;f or uniqueness we must have 
f s ij Vjj ^ dY > 0 

Following Hill (1959), the linear rate constitutive 
equations can be expressed in terms of the Jaumann 
derivative of the Kirchhoff stress, D^. -/Dt, as follows 


Dv 


13 _ 


Ft 


G i3£L 6 kl 


(2.16) 


such that C i;j]a - 0j i]a - G i;jlk = Using the 

relation between s i; j and D T . ./p-t, viz. 

1 J * 


S i3 


^il 

+ a ik v j,k " °ik S 3*k " a 3k e ik 


Dt 


(2.17) 

the uniqueness (stability) criterion takes the form 


Y 


^ ^ G ijkl G io G kl" 2cr i3 e kj e ki + a ij v k,i v k, P dY > 0 


(2.17) 


It is clear that once the instantaneous moduli and 

the stresses o. • are regarded as having assigned values, 

X J 

the criterion (2.17) may be used to determine within what 
limits these values must lie in order to ensure unique- 
ness (stability). 


Finally, like in the classical theory, certain 
restrictions are imposed on the f current ! we 

suppose them to be such that 
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D* 

Dt 


id 


id ~ G ijkl e ij e ki > 0 


( 2 * 18 ) 


for all non-zero stress-rates. This enables one to 
avoid considering the possibility of bifurcation in 
a state of zero stress. *ich (2 . 17) might othermlge 
permit. 
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INELASTIC PLATE BUCKLING 

3.1 INTRODUCTION 

The formulation developed in Chapter 2 is now 
applied to study the bifurcation of thin plates of 
elastic-plastic material. Unlike the case of elastic 
bu-ckling, there have been only a few investigations 
related to the inelastic plate buckling, the most 
notable ones are those of Sewell ( 1963,64) wherein, for 
the first time, the continuum theory of bifurcation was 
applied. It was shown that the bifurcation stress 
(buckling stress) is very sensitive to the local yield 
surface normal; in particular the ’fan’ of possible 
normals at a Tresca corner can give the buckling load 
markedly lower than that obtained using the Mises 
smooth yield surface. Similar observations have also 
been reported by Dubey and Ariaratnam (1969). 

The present investigation is an extension of 
Sewell’s (1964-) work. The bifurcation stress is calculated 
for a simply supported plate in uniaxial compression 
when the stress is at a vertex of the Tresca yield 
surface. The effect of coupled hardening between the 
yield surface facts #iich meet at the vertex is included. 
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A considerable reduction in the bifurcation stress is 
found which, is due to (i) the Tresca comer itself with 
no coupling and (ii) the coupled hardening. 


3.2 FORMULATION QE 1 THE PROBLEM- 


Consider a thin rectangular plate which in the 
current state has edges x-j = 0, a and x 2 = 0, b and is 
bounded by plane faces x^ = + t/2. The plate is simply 
supported on all edges. The current stress distribution 
is uniform uniaxial compression in the x-j - direction i.e.. 


a i3 


(3.1) 


with uniform Q > 0 (Fig. 3.1). Whatever the previous 
stfess and strain history may have been, we suppose 
that the linear rate equations of the following type 
are given at each point.* 


" e n “ 


°11 °12 °13 


^Dt-j ^/Dt " 

e 22 ; 

' = 

°21 c 22 c 23 


D T 22 /Dt 

^ e 33 ' 


O31 c^ 2 c 33 


P T 35 /Dt 


and 2GK, Q 23 = D^/Dt, e 1 3 = 

2G^ 6-j 2 = D 2 /Pt (3. 3) 

There is absent here not only a transverse but 
also in-plane shear interaction with the normal terms. 

It is assumed that the determinant D (say) of coefficients 



IS 



FIG- 3-1 SIMPLY SUPPORTED PLATE UNDER 
COMPRESSION 


in (3.2) is non-zero. Then, (3.2) has a unique inverse of 
the form (2.16) which we write as 


~D^ 1 .j/Dt 
D T 22 /Dt 
_DT 33 /Dt 




C 1 2 
c 22 
C23 




(3.4) 


where a 'bar* denotes the cofactor of the associated 
element in (3.2). For (2,18) to hold, we will have to 
ensure that 


C 11 c 22 



D > 0 


Oil > 0 * 

Gr<l > 0 , 


&2 > 0 


t 


> 0 , 


(3.5) 
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The first of (3.6) is satisfied by actually imposing 
Dtr^/Dt = 0 throughout the plate thickness. This 
condition may be written as 

C 1 3 8 11 + c 23 e 23 + °33 e 33 = 0 (3.7) 

Therefore, the corresponding 'reduced' constitutive 
equations relating the inplane normal stress-rate and 
strain-rates are 



and 







(3.8) 


(3.9) 



c 22 


( 3 . 10 ) 
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(3.5) imply that 




' 11 * c 22 » 


G 3 > 0 ; 


^ *2 

°11 c 22 ~ c * 2 = 


1 


C 11 c 22” c f 2 


> 0 


(3.11) 

The last two of ( 3 . 6 ) can. be taken care of by requiring 
that 8-j 2 ^od 6 2 ^ are small (not necessarily zero) 
compared to 8 ^ and 622 * We now select a class of 
velocity field which satisfy (3.7) and also incorporate 
the usual engineering assumption that the line elements 
perpendicular to the middle surface incipiently remain 
straight and perpendicular to it. This class is 
(Sewell, 1963) 


V' 


U - X 


V - X 


3w 
5 dx^ 

3w 
5 9x/- 


v„ = 


w - x^ ( 


+ \ ( 
2 3 


°1 3 

3u 

J22 

9 v 

'35 

3 x^ 

c 33 

3x 2 

3 

3 w 

°23 

3 2 w 

5 33 

** 'A + 

ax 2 

c 33 

3 x| 


( 3 . 12 ) 


Here u, v, w are arbitrary functions of x.j and x 2 
representing the velocity components of the middle 
surface of the plate. 
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3.3 APPLICATION OF THE UNIQUENESS CRITERION 

The next few steps are similar to those used in 
elastic plate buckling analysis. Here, we essentially 
follow Sewell (1963) and do not reproduce several routine 
calculations. The stress distribution (3»l) and the 
velocity field (3.12) are inserted in the uniqueness/ 
stability criterion (2.17). The effect of transverse 
shear stiffening (i.e„, lj and G- 2 terms) is nearly 
negligible. In the minimisation process of (2.17) it 
can be shown that (Sewell, 1963) the best choice of 
functions u (x., ,x 2 ) and v (x 1 ,x 2 ) is just u = v = 0 
(apart from a possible rigid body rotation). 'Therefore, 
it remains to minimise with respect to transverse velocity 
field w(x-|,x 2 ) of the simply supported middle surface. 
Then, interms of w(x. ,x 2 ), the uniqueness criterion 
becomes 


t 3 a b 8 2 9 2 w d 2 w .2 P 

12 l { [ {0 *1 ( ^5 ) + 2 c t2 py 1^ + °?2 


C 22 


9 w 2 12 3w 2 3 w 

+ 4d (-• ) > - o { — (- — ) - (— ->r) - (■ 

5 3x.j 3 x 2 t 2 9x l 


dx^ 3 x 2 
(3.13) 


) } ] dydx >0 


Satisfying all the boundary conditions, the transverse 
velocity field can be taken in the form 
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m7tX i m%Xn 

w - sin — Sla__ (3.14) 

where m and n are positive integers. Use of (3.H) 
reduces (3.13) to the form 


m 4 * 

t<r> •*, 



12 

r- 

TC~t 


_ (Sr ah- T ~ - (— ) 
2T2 V L a' 12 b ' 


it, 2 n 2 


nt x 2 a 
12 


] > 0 (3.15). 


In addition to ( 3 * 1 1 ) we also require that 

Ci g 'l' 2 ^ 0 


(3.16) 


It follows that, for any given m and a , the left side 
of (3.15) is least when n = 1, i.e,, there is only one 
half-wave perpendicular to the loading direction. 

Also, the thin-plate assumption obviates the need of 

p p 

retaining the terms involving (t/a) and (t/b) in 
(3.15); consequently it is sufficient to express the 


bifurcation stress formula ss 


a = min 


2 2 

fit mb. 2 * 

r [ (-s-) 0 

i2t. 2 a 


+ G-^) 






(3.17) 


It remains to identify the coefficients of the constitutive 
equation (3.2) so that once they are known, the right 
side of (3.17) can be minimised (for a given b/a) with 
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respect to the integer 'm\ However if we relax 4 , the 
constraint that m be integral, then minimisation can 
actually be carried with respect to the continuous ’m’ 
-and the minimum is attained at 


mb 

a 


* 

c 22 
C 1 1 


V4 


( 3 . 18 ) 


Then the corresponding bifurcation stress foimula is • 

[ c* 2 + 2 G 3 + ] (3.19) 


2 + 2 
71 t 


6 b‘ 


3.4 CONSTITUTIVE RELATIONS FOR ELASTIC-PLASTIC MATERIALS 

Here, we refer to Hill’s (1966) constitutive . 
relations for incremental deformation of metal crystals 
by multislip. These relations take into account the 
cross coupling between glide and hardening on different 
si ip -systems. 


The incremental deformation in a crystal is 
treated as the sum of contributions from two independent 
atomic mechanisms on a macroscopic scale Jthese contribu- 
tions are (i) an overall elastic distortion of the lattices 
and (ii) plastic simple shears that do not disturb the 
lattice geometry. The shears have specific planes and 
directions each associated pair of which defines a 

+ can be justified in many cases as shown in Sec. 3.5. 
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glide system'. -At each stage of a loading programme 
the particular 'active' systems are considered to be 
uniquely determined by the current stress-distribution; 
its shearing component in a direction of possible glide 
must attain a critical positive value that depends on 
the history of overall straining and on the particular 
system. luring incremental deformation each system 
hardens at a definite rate which may be a function of 
every rate of shear, typically y say. 

In an elastic element we take the rate equations 

to be 

ij K ijkl — (3.20) 

where the moduli have the symmetry of in (2.16): 

K ijkl = K 5ikl = K i 3 -ik = K klij (3.2l) 

Ihe plastic part of the overall strain rate 6- - is 

^ d 

expressed as 

ij ~ Y a m ija (3.22) 

Greek suffixes will range over 1 to IT, for some specified 
IT, with summation implied when they are repeated unless 
otherwise mentioned. Each y a is interpretable as a 

shear rate in a direction with unit vector components t- 

■ xoc 

on a glide plane with unit normal components n ifr , with 

J'VV 
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m ij a 2 n 3 + n i^a (3.23) 

Thus nijL j a m i; j a = \ (no sum on a) ( 3 . 24 ) 

Accordingly to Sewell ( 1 97 1 ) the precise single crystal 
interpretations do not have to be taken over directly. 
Instead, the connection with more classical plasticity 
can associate the N tensors f2 with the unit 

normals to the facets of a singular yield surface having 
a vertex at the given stress point. We shall illustrate 
h = 2 with the Tresca criterion. When U = 1, equations 
appropriate to a smooth yield surface can be recovered. 
Following Hill (1966) and Sewell (1971), the right hand 
side of (3.22) is expressed in the form 

-1 D 

Y a m i;ja = m klp' m ija 

in which the h^ are the elements of a given N x 3ST 
'hardening' matrix such that 

positive definite (3.25) 

Finally the constitutive relations for an elastic 
plastic element become 

-1 DT kl 

6 ij = ^ijkl + m iD« m m ) 5t (3.26) 

This is the form of the constitutive law from which we 
can expect to identify coefficients in (3.2) and (3.3). 
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We adopt the Tresca yield criterion which»in the 
plane a ^ - 0 of principal stress space, is shown in 
Fig. 3.2. We assume that the principal axes of and 

m ij2 coin cide with those of stress. Then it follows that ; 



FIG’ 3*2 TRESCA YIELD SURFACE 



We further assume that the elastic part of the material 
response is isotropic with Young's Modulus E and Poisson's 
ratio v. Then \i = E/2 (1+ v) is the positive elastic 
shear* modulus. Further, the 2x2 assigned hardening 
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matrix is written as 

r i n 



in which. 

h > o , - i < n < i 


(3.28) 


(3.29) 


in order to satisfy (3.26). The parameter * represents 
the amount of ’coupled hardening*. It may be mentioned 
here that in the investigation by Sewell (1964), the 
hardening parameter h is twice the present ’h 1 when 
N = 1 and the yield surface normal is So far as 

is known, only three types of hardening matrix have even 
been proposed, and all are special cases of (3.28): 


(i) 9 = 1 f isotropic hardening (Taylor) 

(ii) "H = 0, independent hardening (Eoiter) 

(iii) 9 = 1/2, kinematic hardening (Budiansky and Y/u). 


The relations (3.27) and (3.28) and the isotropic 
form of are now substituted into (3.25). The 

constitutive equations then take the following form 


e 11 


e 22 

1 

" E 

e 33 

1 »- 


0 


1 + 


— v— 


— v. 


u v. 


1 + 7 } 

6 

2(1+*) 

1 

2(1+ 1 


1 + 


2 ( 1 +*) 

0 

2(i-* 2 ) 

TJ0 


— V — 


— V — 


2(1+.*) 

TfQ 

2(1-t> 2 ) 

0 


— V 


1 + 


2(l-n 2 ) 2(l-rj ) 


Dt 


11 


Dt 


Dt 

Dt 


22 


D 

Ut 

T3.30) 
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and 2[x 


- 2n S 13 = 


’ 2 ^ ®1 2 = 


where 0 — E/2h. From (3.30) are oan easily calculate 
the starred moduli using (3.10): 


(l + 


2(l--o ) 


„ e e 

C = ( V + ) 

12 D 2(i +ri ) 


(3.31) 


C 22 D ^ 1 + 


1 + T) 


S * ... E 

&„ = > c . i C — c* = 

3 2( 1 + v ) 11 22 12 K 


0 e 

where D = (l + -) (l + 

(l+T)) 2(1-- 


■) - (■ 


‘rV! -TJ 2(1+tj) 

All the relations (3.1 1) on the moduli are therefore 


satisfied. We also need 


; l2 + 2 & 3 = 5“ ( v + 


e e 

— ) + — 

2 ( 1 -HT).) 1 + V 


(3.32) 


and it is again -verified that (3.10) is satisfied, too . 


3.5 CRITICAL STRESS AT BIFURCATION 


The bifurcation stress formula appropriate to a 
comer on the Tresca. yield surface is obtained by inserting 
(3.31) and (3.32) in (3.17); it can be written as 




29 


12 ab 2 1 . v 2 a 

~ % £ t 2 ^ ~ ^ I'rp (m,e,n~ » v ) (3.33) 


where the function of the indicated variables is 





2 

) 


2-2rj 2 + e 

2( 1— T) 2 ) 


2 v + 2 -nvt- 9 

+ 2 { 

2(1+ T] ) 


“h 


D 

1 + v 


} 



1 +T}+0 
+71 


(3.34) 


in which D itself is a function of (9, v , rj ) as listed in 
(3.31). The formula appropriate to a smooth yield surface 
N = 1, essentially corresponding to the case of Mises 
yield surface, has been given by Sewell ( 1964 ); it can 
expressed as 


2 2 

12 ab 1- v a 

) = Min ( m , 0, v ) (3.35) 

r t E D 


where the f unction F^ is 


M 


mb 2 Q a 2 20 

("a - ) 0 + 5 ) + ( m5 ) ( 1+ “ )+ 2 + 

5-4 v 9 

[ 1+ (— ) 9 ]/ (1- v 2 ) 


9(l-2v) 

3(l+v) 


(3.36) 


The aim is now to find out whether the new formula (3.33), 
appropriate to a yield surface corner with hardening 
coupled according to the factor t) , can predict a 
different bifurcation stress from that given by formula 
(3.35) corresponding to the Mises case. A numerical 
investigation has been carried out. Throughout we let 
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the value of the Poissions ratio v= 1/3. The right 
side of (3.33) and (3.35) were plotted in the range 
® i ® ^ for ten aspect ratios in the range 

a/b = 0.5, 1, 1.5 - - - 4.5, 5. It was found adequate 
in (3.33) to fix on five values of p within ( 3 . 29 ), 
namely p = -0.9, -0.6, 0, 0.5 and 0.9. This makes six 
plots in all for each aspect ratio. Pigs. 3.3 - 3.5 
show such plots for a/b = 1.0, 1.5, 3.0. Also P T and 
Pj^ were plotted as functions of a/b for a given value 
of 9 and p (in case of P^) . Pig. (3.6-3. 8) show some 
sample plots, which are the results of minimisation of 
Pj and with respect to positive integral values of m. 
The m value for an aspect ratio seems to depend on both 

8 and p . Por example, for a/b = 1.2 the m is 1 for 

9 = 0 (elastic case), it equals 2 when 9 = 10 for the 
Mises case as well as for the Tresca case (p =0.9) but 
takes on the value 1 for the latter case if p = -0.9. 


Next we minimise (3.33) and (3.35) with respect 
to continuous m. The associated formula for the Tresca 
case is 


2 2 

12a b 1 - v * 

TTr < — > ( 9 ^ .*> 


71 V 


, T,* 

when -P ri = 


2(1— v^) ^ (2-2p 2 +9)(l + pf9) 

5 2(1 - P 2 )Vl + -p) 


(3.37) 

V 2 2v +2 p v+ e D n 

} + . .. . - - — » - ■ + J 

2 ( 1+P') 1+v 

(3.38) 
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The minimisation can actually be carried out ■with respect 
to the parameter mb/a and as in (3.18), the minimum is 
attained at 


mb 

a 


^ 1 + 0/(1+ n) 1/4 

1 + 9 / 2 ( 1 .^^) ? 


(3.39) 


The corresponding formula for the Mises case is 


12 cr b 2 1 - v 2 


~~T'T 

TC t 


( 


E 


) ~ % (0, v ) 


(3.40) 


where 


r . 29 x 9 1/2 0(7-2v) 

[2 (1+ r~ )(l+ -~) +2 + ~r N ] 

f 5 - T_— — — hit.'-) - (3.4i) 

5-4 v ■ 0 

[ i + ( — )e] /(i-v 2 ) 


The analytical minimum is obtained, in this case, at 
the value 
mb 


6 + 46 1//4 

= [- J - 


6 + 0 
* 


(3.42) 


The functions Ejj and were plotted for many possible 
combinations of the variables 0 and p . The horizontal 
lines in Fig. 3. 6-3. 8 are, in fact, these functions. 

When compared with the graphs of integer minima, it was 
found that there was not much difference between the use 
of integral or continuous m in the minimisation, at least 
for a/b > 3. This might have been expected from what is 
known about the classical long plate approximation in 
elasticity, but it was necessary to check in the present 



context also. It is, therefore, adequate to concentrate 
attention upon six curves (five functions F* and the 
function Fjj ) shown in Fig. 3.9 obtained by minimising 
the right side of (5.37) and (3.4O) with respect to 
continuous m . 

It is to be noted from Fig. 3.9 that all the 

curves for Tresca case lie below the Mises curve as 8 

increases from zero. Infact, an increase In 9 from zero 

represents the fall -off in tangent-modulus from E, as 

will be shown later. Even for independent hardening 

(r) = 0) for the two yield surface facets, the function 

Fji drops considerably below the Mises curve (F^) and as 

■H increases into positive values the gap widens at an 

increasing rate. Asp-*- 1 , it is clear that a dramatic 

reduction in the ratio F(j» / F^ takes place even while 9 

is very small. We cannot actually put t) • = 1 as it 

will violate (3.29), but for h = 0.9 We find that 

V F M = 1/1.6 is achieved by the time 8 has reached 

about 3.0 . In the absence of clear information from 

physics about the value of P , we have investigated the 

p 

effect of all "H < 1 . It turns out, as seen In Fig. 3.9 

* 

that the above mentioned sensitivity of F^ to variations 
in positive value of T) does not occur for negative tj . 
More-over, the F* curves for > 0 do not intersect 
each other, in contrast to those for 9 < 0 . 
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The minimising' "values of mb/a, given by (3.39) and 
(3«42), are shown inFig.3.10, as a function of 9 , There 
is something special when the coupled hardening has a 
value p = 0.5; the matrix of moduli in ( 3 . 30 ) becomes 
isotropic instead of anisotropic, just as for an elastic 
solid. This is why the horizontal line at mb/a = 1 is 
obtained f or p = 0.5 in Fig. 3 . 10. It implies the 
continuation into the plastic regime of the elastic 
property that for plates of integral ratio a/b, the 
buckling half -wavelength and the aspect ratio coincide 
and are independent of the amount of strain (represented 
by 9 here) before bifurcation. Fig. 3.10 also shows 
that the critical value of mb/a falls with increasing 9 
when 0.5 <p < 1 in constrast to its rising behaviour 

for - 1 < p < 0,5 and in the Mises case. 

3,6 NUMERICAL CALCULATION OF BIFURCATION STRESS 

It has been shown in the last section that sub- 
stantial differences do occur between the right sides 
of the bifurcation stress formulae (3.33), (3.35) and 
their respective lower bounds. As shown in Fig. 3.9, 
the error involved is negligible if a/b > 3 . We now 
wish to compare the bifurcation stresses actually 
predicted in this range of a/b by the formulae (3.37). and 
( 3 . 40 ), in situations corresponding to an actual test. 
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In order that any F^, , 0 point or Fj/[ , 8 point in. 

Fig. 3.9 can actually represent a bifurcation value of 0 
(or of the hardening parameter h if E is constant), the 
strain history must have been such that the corresponding 
bifurcation stress a required by (3*37) or (3.40) be at 
a point of the stress-strain curve which is compactable 
with this value of 9 (or of h ). 

To achieve this, and so to obtain and compare 
some numerical values of the bifurcation stress, we make 
some assumptions which will enable us to determine the 
left side of (3.37) and (3.40) as an increasing function 
of 9. Then these curves can be superimposed on the 
decreasing curves of Fig. 3.9; the corresponding intersections 
will determine actual bifurcation values of a and 9. 


Suppose that the instantaneous tangent modulus 3^ 
in compression is adequately expressed as l/c-j-j from 
(3.8), by inserting Di^/Dt = 0 there to give 6^=0^ Dtj^/at. 
From (3.10) and (3.31) it follows that 


E ® 

“ 1 + n 


(3.43) 


for the Tresca case. For the Mises case, it is' given by 
( Sewell, 1964) 




(3.44) 
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It is to be noted that when p = 0.5 , the two tangent moduli 
are the same functions of 9. This is the same amount of 
coupled hardening which results in some simplifications 
mentioned in connection with Fig. 3.10. 

It is clear from (3.43) or (3.44) that as progress 
is made from yield along a stress-strain curve with 
monotonic ally decreasing E-f., the parameter 9 will mono- 
tonically increases from zero. We neglect any variation 
with 9 which the width /thickness ratio (b/t) may have. 

Then, with any given stress-strain curve (and fixed v ), 
the left side of ( 3 « 37 ) and ( 3 . 40 ) is expressible as 
a decreasing function of B^.. For each fixed p it will 
be an increasing function of 9 , from ( 3 . 43 ). For 
example, consider a two-segment stress-strain curve 
represented by the Ramberg-Osgood relation of the type 




+ a 



(3.45) 


where a, and n are material parameters, ey is the yield 
strain which equals cr /E, a is the yield stress and e 
is the compressive strain such that de/dtf = l/E^. Hence 



a n -1 

= 1 + an ( -*r~) 


and therefore from (3*43) we get 




(3.46) 


(3.47) 
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Now, the left side of (3.37} becomes 


12 ob 

“2 2 
TU “G 


in which 




an(l+rj) 
12(1 -^) 


1/(n-l) 

] a* 


(3.48) 

(3.49) 


S'or the purpose of numerical calculation we take a =0.1, 
n = 8 (Needleman, 1973) and consider variations in c*wb.ich 
is identical to varying the ratio of o> /E or the ratio 
b/t (or both) because v is fixed. For the sake of 
illustration, three values of a* (i.e. 2.5, 3.0, 3.5) + are 
taken. For each o* the rising curves (3.48) are obtained 
for five values of T) = - 0.9, -0.6, 0, +0.5, +0,9 and 
are plotted in Fig. (3.11-3.13). We recall from (3.44-) 
that the curve (3.48) for p = 0.5 will also be for the 
Mises case. The curves from Fig. 3.9 are repeated in 
each of the Fig. 3.11-3.13; the corresponding intersections 
are bifurcation points. The common values of the two sides 
of (3.37) and (3.40) at bifurcation, as determined from 
these intersection curves, are shown in Table 3*t. 


The conclusion to be drawn is that the occuranee 
of a corner on the yield surface at the local stress point 


+ go get an idea about the plate thickness, we note that 

for steel (c^.= 2500 kg/cm , B=2.1 x 10° kg/cni ) as a* varies 
from 2.5 to 3.5, b/t varies from 44 to 52. 
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can cause a lowering of the bifurcation stress when the 
coupled hardening coefficient does not assume a negative 
value. For r)2 0, the reduction does exist. In particular, 
the lowering can be present even when hardening of the 
facets making up of the yield surface corner is uncoupled 
( T) = 0), and a greater reduction in the bifurcation stress 
can occur if the coupling extends much beyond the value 
h = 0.5. It is to be noted that this reduction has been 
achieved even though the shear modulus retained its elastic 
value throught. 
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FI&-3-12 BIFURCATION STRESSES CORRESPONDING- TO STRESS-STRAIN CURVE (EQ.3-45) 
FOR « = CM, n = 8, crt3*0 






CHAPTER 4 


ELASTIC PLATE BUCKLING 


4. 1 INTRODUCTION 

The study of the buckling phenomenon in thin 
elastic rectangular plates has been carried out by 
many investigators in the light of small deformation 
theory • The aim of this chapter is to reexamine this 
problem under the present framework of uniqueness and 
stability, outlined in Chapter 2* The plate material 
is assumed isotropic in the ground state but is allowed 
to develope anisotropy during deformation; the material 
becomes orthotropic in the current state the stability 
of which is being investigated. The plate is assumed 
simply supported on all edges. 

The method of analysis is similar to that in 
Chapter 3; the difference lies in the constitutive 
equations used. The plate is supposed to be made of 
a general incompressible elastic material by adopting 
the strain-energy function of Ogden (197 2). The energy 
function contains a parameter a by varying which a large 

..V '■ V D ; \ y ' 7 u V: 

variety of solids can be generated* Incidently, a = 2 
will correspond to a neo— Hooke an material* One of the 
advantages of using such functions is that the results 
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(i.e. the critical stress) can be expressed inte rm s of 
moduli (namely the shear modulus in the present context) 
in the undeformed configuration. The rate constitutive 
behaviour of the material is taken in the form proposed 
by Hill (1970) which is available in a form simpler and 
more convenient than any other form involving strain- 
invariants. 

4.2 CONSTITUTIVE RELATIONS 

The principal Cauchy stresses Oj , a 2 * may 
regarded as functions of the three principal stretches 
a.j , a 2 >a^ measured from a ground state with respect to 
which the material is isotropic. The principal axes of 
the stresses are coincident with the axes of the Eolerian 
strain ellipsoid. The principal Kirchhoff stresses are 
defined simply by 

•'l = J °1 - x 2 - 3 « 2 . ^ = j 

where J = a-ja^. If the solid is hyperelastic with 
strain energy W per unit reference voltane, we have 


3W 


dW 


dW 


X 1 = a 1 


9a 1 ’ ^ 9a 2 


a 3 3a, 


(4.2) 


The rate of change of principal Kirchhoff stresses are 
given as 
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3Xj 


~ ®1 1 + ®22 + a 3 ~ e 33*. etC ' 
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cdMl 


UBgAgy 


A fififtl Q 
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where a,^^- S 11 t a 2 ./a 2 = e 22 ,i 3 /a 5 = e 


'33 


are "the normal components of "the Euler ian strain— rate 
®ij 0n °f "the Eulerian ellipsoid. Then, the 

incremental relations in terms of the Jaumann derivative 
can be expressed as (Hill, 1970) 


BT 1 1 l , 9 a 

3>T 8a, ' l ' a 1 a 1- ,e i1 + 8aJ (a 2°P e 22 + 8^ C^o, Je^.eto. , 


Pt 12 = 
Dt 


2 2 

a,, + a 2 

a 2 _ a 2 (° 1 " a 2 ^ e 12 > *1 ^ a 2 
1 2 »etc. ( 4 . 4 ) 

a 

a 1 daj’^ Cr i” a 2^ e l2 » ( a 1 = a 2^ 


which is in a form similar to (2,16) ideally suited for 
using in the uniqueness/stability criterion. The instan- 
taneous moduli are then just the coefficients of the 
components e^* in ( 4 . 4 ), 

In case the solid is isotropic in the reference 

configuration, the potential W is a function of special 

choice of strain invariants (e.g. Green and Zerna, 1968) 

or of principal stretches a 1 as considered by Ogden (1972). 

We consider the following strain-energy function for an ' 

incompressible elastic material 
a a a , 

a., + a Q + a, - 3 .. . x 

W = L (~ 2 '- - 2 ) (4.5) 

r . a ; ; /■■■-. : ; ,v 

where oc is a real number and f* the mechanical parameter 


such that 
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11 “ > 0 (4.6) 

is ills shear modulus in the undeformed configuration. 
The principal Cauchy stresses for an incompressible 
solid are expressed as 

1 -j 9W . 1 9W 

a^a*2 ^ 2 a^a^ 9a£ 5 ^ a 3 

(4.7) 

where p is an undetermined quantity, For the function W 
defined by (4.5) , the relations (4.7) become 


°ij - Pa^ = “ ^ a 2 = ~ ^3 = P (4.8) 

such that 

a 1 a 2 a 3 = J = 1 (4.9) 

which is the condition of incompressibility. Since the 


current state of the stress is the uniaxial compression 
along the x^ -direction (i.e., = o^ = 0) , we have 

from (4.9) 



= aT 1 / 2 


( 4 . 10 ) 


so that we get 



a -a/2 
3*1 — 3 ^ 



(4.11) 
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4.3 APPLICATION OP UNIQUENESS CRITERION 

The first term in the uni quenes s/stability 
criterion (2.17) can be expanded in the following form 
in view of (4.4) 


BT ij p 1 a aw 2 

Dt id ~ C ijRl ®io e ih = “al, (a 1 9^ 6 i 1 1 

ia aw o 1 a 


a 3 a -j 3 b. 2 ^ ^ a 1 a 2 ^ ^ a 3 ” 

2 A 2 d 2 w 

a 1 9 a 2 3a 5 22 53 a g aa 3 3a 1 53 11 
2 a 2 w 2(a|+a|) 2 

+ ^yiy eii ® 22+ uf^r 1 Co 2-^ e 25 


( a 3 + ai ) , s p 2(a^+a^) ^ ^ P 

T“5 5a ( o; 2 - 0 i ) 6^ 1 + —75 — 7T ( o’ . — ox ) 8., 


l 3 " ®1 


) v 3-1-31 -7^7 


J 1 u 2' °12 


( 4 . 12 ) 


which after using (4.5) becomes 


BT iD 


a „2 


= *“ ^ efi + »*“ 4 e | 2 + 2 ( g - ) ( 0 l - cr 2 ) e * 2 


2 2 
3^ 3«2 


(af + a| ) , 


a 1 ~ a 2 


(a?- A 


-V e; 


(4.13) 


where a^a 2 a^= 1 and o-pc^o^ are given by (4.8). Also 
other terms in (2.17) get simplified to 


/ 2 2 2 2 2 

— ®^j( 8^ + S-| 2 + ®1 3 *” “12 “ « 1 j + 28 -j 2 ^1 2 2 ®1 3 ®-| 3^ 

( 4 . 14 ) 
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so that (2.17) can now be written as 

DT. 


/[ ™. 
Y I>t 


13 


e i3 -°1 < a 


11 


e 2 + p 2 „ 2 
fc 12 fc 1 3 ~ *12 


4) 


+ 2e 1 3 *i 3^ ] dv > 0 
where the first term is ( 4 . 13 ) itself. 


13 + 26 1 2 W 1 2 
(4.15) 


In the current state, the plate has dimensions 
a x b x t. 'The plate is simply supported on edges 
X 1 = °> a and x 2 = 0 , b while the faces x^ = + t /2 are 
free from applied traction - rates i.e., 


• « 

s 33 ~ s 31 = s 32 =0 on x 5 = + t/2 ( 4 . 16) 

which are same as ( 3 . 6 ) . The boundary conditions can 
be taken care of in a manner similar to one in Chapter 3. 
The velocity field, infact is also identical to ( 3 . 12 ). 

In order to avoid the repetition of some unnecessary 
steps, we consider the following subset of the velocity 
field. 




3w 




V 1 = 

u 

" x 3 "a^ - u 

1 

M 

< 



II 

C\J 

f> 

u 

dw 

- X, = u 

3 3x 2 ~ 

- x 3 

(4. 

17) 

<} 

* 

II 

w 






w 

are the middle 

.surface displacements 

and 


are functions of x 1 , x 2 only. Again in the minimisation 
process of (4.15), it turns out that the best choice 
of the velocity functions u and v is u = v = 0 . 
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Then the condition (4.15), for buckling not to set in, 
can be very much simplified to yield 

?2 { aa 1 \ 1 + “ & r /2 ^22 + 2(a 1 “ a^/ 2 ) ( ) wf 2 

a 1" 1 

f cc —(X / 2 1 2 p 2 2 

+ (a,- w n - w#t1 - w 12 ) ]dx-j dx 2 > 0 ( 4 . 18) 

Since the plate is assumed simply supported on all edges, 
we take, as in (3.14) 


mitx-j mmxp 

w = sin sin 


-v (4.19) 

a D 

when m and n are integers. The use of ( 4 . 19) reduces 
the condition ( 4 . 18) to requiring 


m 


4 [««?- (a? - af /2) ] + ( \ ,* + ( S) 2 m 2 n 2 


<4^ < 

a(-1 


2 2 

*1“*! \ , m a / a -a/2. s „ , 0 . 

” — ) + 3 “2“^ (a 1 -a 1 7 ) > 0 (4.20) 

TC w 


Por a given a/b, t/a and a , the critical stretch a^ is 
obtained by minimising (4.20) with respect to integers m 
and n . As expected the minimising value of n is always 1 ; 
the modified form of ( 4 . 20 ) is then 

»+[ - (4-aT“ /2 )] + %-) 4 “aT“ /2 + m 2 ( 


2 2 
a m 


+3 c -^2 (a ? - 4“ /2 5 •! > 


0 


(4.21) 
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4. 4 DISCUSSIONS 

Once the plate material is specified (through 
parameters p and a in (4*3) ) and the plate dimensions are 
given (the aspect ratio a/b and the ratio t/a), the 
expression (4. 21) can be minimised with respect to m to 
yield the critical value of the stretch a-j at bifurcation, 
.Having known the critical stretch, the critical stress 
can be calculated using(4. 11 ). It may be mentioned here 
that the negative values of a are also admitted because 
p can be negative too to make the shear modulus Pa 
positive. 

Six values of a were taken: a = -2.0, -1,5, -0,5, 
0,5, 1.5, 2,0, The plate thickness/width ratio was 
maintained as 1/30 or 1/100. Dor each set, the aspect 
ratio was varied from 0.4 to 5.0 and in each case (4, 21) 
was minimised with respect to the integral m to yield the 
critical stretch (and hence stress) , Dig. 4.1 and Dig, 4. 2 
represent the variation in the buckling stress with 
aspect ratio for t/a = 1/30, 1/lQO respectively, when 
a = 2.0 (i,e,, nea-Hookean material). Dor other values 
of a, it is seen that the nondimensionalised critical 
stress increases as a decreases from 2.0 (excluding the 
value a = 1.0). The increase is, however marginal 
for t/a = 1/30 as can be seen from Table 4.1; and is 
almost negligible for t/a = 1/100 • 
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TAELS 4.1 

VARIATION OB' CRITICAL STRESS WITH a EOR t/a = 1/30 


a/b 

a 

Critical Stress 

a /pa x 10" 

6 

2.0 

1.5 

0.5 

-1.5 

0.4 


7680 

7696 

7727 

7791 

» 

1.0 


3660 

3661 

3663 

■3668 

1.6 


3844 

3846 

3850 

3859 

2.4 


3782 

3785 

3786 

3789 

5.6 


3700 

3701 

3702 

3711 

5.0 


3600 

3661 

3663 

3668 


It is interesting to note that though a large 
number of solids can be generated by varying a, tlie 
critical stress is not really sensitive to the variation 
in a . This is certainly true for the range of t/a 
considered here. No conclusion can be drawn for thicker 
plates because the analysis itself would require 
modification. 
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FIG- 4 1 VARIATION IN BUCKLING STRESS WITH ASPECT RATIO (t/a 
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CHAPTER 5 


CONCLUSIONS 


This thesis has been concerned with the investigation 
of the buckling phenomenon in thin elastic-plastic and 
elastic plates under finite axial compression. An attempt 
has been made to explain the so-called paradox in plastic 
plate buckling theory; the paradox is that the experimental 
buckling loads have been markedly lower than the theoretical 
bifurcation load obtained from the Prandtl-Reuss theory 
using the Mises smooth yield surface. 

The bifurcation (buckling) stress has been calculated 
for a simply supported elastic -plastic plate, in the case 
when the stress is at a vertex of the Tresca yield surface. 
The effect of the coupled hardening between yield surface 
facets which meet at the vertex has been included. Pig. 5.9 
shows how the bifurcation stress depends on the tangent 

modulus of the compressive stress-strain curve. the 

parameter 9 is the reciprocal of the ordinary hardening 
normalized by Young's modulus E, so that as 9 increases 
from zero the tangent modulus decreases from E. Pig. 5*9 
also shows how the critical stress depends on the measure 
U of coupled hardening. This dependence has been found 
to become more sensitive as tj increases from 0 to 1 • 
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Pig. 3.11 to Pig. 3.13 show how the Tresca bifurcation 
stress falls below that of Mises for some particular 
value of various parameter in the stress-strain curve 
which is of the Ramberg-Osgood type. It has been 
concluded that the reductions of order 10 to 30 percent 
may not be unrealistic when Mises is replaced by Tresca. 
Two separate new effects, infact, contribute to the 
reduction obtained. One is the comer itself when there 
is no coupling (r) = 0) ; the other is coupling. Much of 
the previous literature is concerned with theories which 
embody a shear modulus which drops below its elastic 
value. The present investigation, however, has not taken 
this feature into account. The shear modulus does appear 
in the bifurcation stress formulae, but it is not the 
only parameter there, and the stated reduction have been 
obtained by using the above two effects together with 
an elastic shear modulus. 

Por the inelastic plate, simple probabilistic 
formulation has been given to predict the buckling stress 
when the random variables follow; the normal distribution. 
Por a desired reliability (or the probability of failure) 
and an assumed variation in material and geometric 
parameters, the buckling stress has been calculated and 
compared with those obtained in Chapter 3; the results 
based on the reliability analysis are always on lower 


side 
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All the previous investigations concerning 
elastic plate buckling have presented results in terms 
of 'current 1 moduli of the material; these moduli, in 
turn, are related to the moduli in the undeformed state 
of the body and these relations are obtained by solving 
the problem of finite deformation from the initial state 
to the ground state. This has been obviated in the 
present investigation by considering the Ogden's strain- 
energy function; the critical stress, thus obtained is 
in terms of ground state shear modulus. The material 
is initially isotropic but is allowed to develop 
anisotropy through deformation. Although, only incompre- 
ssible materials are considered, the analysis can easily 
be extended to utilise the strain energy function for 
compressible materials. 
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(a) m=1 



FIG- 5-1 BUCKLED FORM OF PLATE SUBJECTED TO 
IN-PLANE COMPRESSIVE LOADING 
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APPENDIX A 

A RELIABILITY APPROACH TO INELASTIC BUCKLING 


Consider a structure (or structural component) 
whose resistance is a random variable X , subjected to 
a random load Y with known density functions of (x) 
a^d fy(y), respectively. The failure is then defined 
as the occurance of the event ( X < l) and the associated 
probability of failure is (Fig. A-l). 


P f = P ( X< Y) 

°° I 

^ [ f f (x) dx ] f (y) dy 

O O * 

oo 

= / F (y) f (y) dy (A. l) 

Q & X 

in which Fy( ) is the cumulative distribution function 
of Z. When the load is I = y the probability of 
failure is F-^y). If x and y are normal variates with 
mearsX and f and standard deviations S.j and Sy, the 
failure may be defined as ( X-Y < 0) and the calculation 
of Py is simplified considerably. In this case (X-Y) 
is also a normal variate with mean ( X - Y ) and variance 

p p 

( + S^). Then the probability of failure is given by 

P f = P (X - Y < 0 ) 
and the equation (A.l) becomes 
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X - t 

(sf + s|)V2 ] 


in which 0(z) is the tabulated cumulative probability 
of the standard normal variate (Z) for Z <. z, ■where z is 
given by 


z 


X - Y 

-~ 1/2 

( s x + ^y 


(A. 2) 


In (A. 2), z defines the quantum of probability of 
failure. Once the probability of failure is known, the 
reliability is given by 

P r = 1 - P f 

Now we apply the reliability theory to the 


problem of inelastic plate buckling in Chapter 3. We 
have the strength-stress relation for the inelastic 
plate given by (3.33) as 


1 - v mb 2 

( >[ ( 22 ) 

D a 


,2-2 p 2 +e, 2 +2 71 V+ 9 D 

( 3-) + 2 ( " ' » '■ -f — 

2(i-if) 2(1+ n) i+ v 



1 i+n 


) 3 > 


12b 2 1_ v 2 

( T } 


From the above equation we observe that the expression 
on the left, henceforth called by K c , denotes the 
resistance of the plate that has to be overcome If the 
buckling has to set in. The parameters 6, b, a have a 
random variation with means 9» a and the variations 
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V e’ Y b’ Y a res Pectively. 


We denote the expression on the right as K g , since 
it is the actual stress that will develope at the time of 
buckling and is given by 

2 p 1 

• i2tr i-\t fl & iox2 r . 

K s = a -5- (— ) = [ — _! — >-i hLihl 5 

* B an(l+r})- E ^ t 2 

r e I. . 

an(l+ p) 

involving a , E as random variables, let c7y, B be 


the mean values Y cr y , V B be the variations, respectively. 
So the reliability against buckling is given as 


z 



f s 2 

v s kc 


U.3) 


where K c is the mean expected strengths and K g is the 


mean developed stress at buckling, and S^ c , are the 
respective standard deviations. For a given value of 
Pf the value of z can be obtained from the standard 
probability tables (Appendix B) and vice versa. 


Since each of the variables in K c and K g is 
associated with a statistical distribution (assumed 
normal here), the over all distributions of K c and K g 
can be obtained by knowing the coefficients of variations 
of each of the variables. An approximate expression for 
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standard deviation (and hence variation) is given by 

2 21 dK o , 

s ko “ ? . > x 2 


s te 

= E 
i=1 

~ 

s 2 

m 
= y 

0K S 

s ks 

i=1 

dX i 

where K c 

= f (x 

1 > x 2 


&1’ 


Using the above formulae the following relations can be 
derived? 

ks ° (u-1) 2 8 


where 


V 2 * 
a 


+ vf + 4 - v £ + V 2 


2 ? 2 
v e = 4 +v h 


S + 4V h + v ay 

.2 


Also 


1-v 2 2 2 2^2 2 +e 1 2 2 + r)+ 9 

( S; ) { t » fe) ) + ( T> ( -TT7- ) 


■'to = ( s : 5 { t * 


2 v 4-2n v+e 2 p 2 1 ^ 

+ (__ ) ] V D + ( 2(i.x" 2 ) ) + 1+ U 

*1 O ^ 2 ^ 


l —] 2 e 2 v 2 + tx 2 (!±^ e > - 2 x 2 (l^i -)] 2 v 2 > 

( 1 + u ) ' 0 - 1 - t »2 1 + U - X 


x^( 1+t}) - 0 - i-n2 i+ n 

mb 

where we have taken x = ~ . As already stated for a 

min im um with respect to continuous m, we have from 
Chapter 3, 

2 (i-r)) (i+ u+e ) . 1/4 

x = [ ' 3 

2 

2-2n + e 
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cn +* + u2 0 r (W) (l-2rf) 

so that — - [ 

4X (2-2ti 2 +0) 2 *5 2(1 -h) (i-frj+e) 1,5 ' 
Similarly tne mean values K c and K c can also be defined 
taking the mean value of each variable as 


= [ 


an(l + rj ) 


•] 


Having known V kc and V ks , a criterion to find the 
buckling stress can be obtained from the equation 
(A. 3) as 


K s T ks> -2 KsK 0 S v£ 0 E* = 0 (A. 4) 

The above equation is a function of 0 only, that will 
determine the stress at buckling (refer Figs. 3.11-3.13 
for intersection curves at various values of 0). For 
convenience, we have taken a constant coefficient of 
variation (V) for all random variables. So solving 
(A. 4-) for the critical value 0 Q at which buckling occurs. 


2 „2 


the buckling stress can be obtained as 

k! = 


an(l+r)) 


Fig. A- 2 to A-7 give the variation in the buckling stress 
with o . Here we point out that a is a parameter which 
represents the material characteristics (cry, E) and 
also the thickness of the plate. For z = 3 (Ff = .00135) 
and for z = 5 (P^=0. 000000 3), the buckling stress can 
be readily found, for variations of 0.01 and 0.05, 0.10. 
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The above values are compared with the corresponding 
values in Chapter 3 for 2.5, 3 and 3.5 (Table A-l). 

The equation (A. 4 ) can also be used with 
a dvantage to probe the sensitivity of the reliability 
to the plate thickness (and hence on the cost of plate) 
for a given overall coefficient of variation and the 
hardening parameter 0, as shown in Figs. A-8 and A-9. 
Such curves are helpful to a designer, so that he can 
select a plate thickness, depending on the reliability 
sought. 



71 


TABLE A-1 

COMPARISON 0E THE CONVENTIONAL BIPURCAIION STRESS 



WITH PR03A3LISTIC 

STRESS FOR 

^=0.05, 

z=5 




=2.5 

■if 


* 



a 

=5.0 


=3.5 



Conv, 

. Probl. 

Conv, 

Probl. 

Conv, 

Probl. 

-0.9 

3.6 

2.15 

3.77 

2,20 

3.90 

2.22 

-0.6 

3.3 

2.10 

3.88 

2.17 

3.75 

2.20 

0.0 

3.1 

2.05 

3.40 

2.15 

3.58 

2.16 

0.5 

2.88 

2.00 

3.20 

2.10 

3.45 

2.11 

0.9 

2.63 

1.85 

2.90 

2.00 

3.15 

2.08 


Conv.. : Conventional method 


Probl. : 


Pro babli stic method 






0 1 V5 2 2-5 3 3-5 

cr 

FIG-A-7 VARIATION IN BUCKLING STRESS WITH <r* 
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APPENDIX B 


STANDARD NORMAL PROBABILITY TABLE + 


2 

f (z) 

Reliability 

R(z) 

Probability of 
f ailure 'Pf * 

0.0 

0. 398942 

0.500000 

0.500000 

0.1 

0.396952 

0.539827 

0.460173 

0.3 

0.381388 

0.617911 

0. 382089 

0.4 

0.368270 

0.655422 

0.344578 

0.6 

0.333225 

0.725747 

0.274253 

0.7 

0 . 31 2254 

0.758036 

0.241964 

0.8 

0.289692 

0.788145 

0.211855 

0.9 

0.266085 

0.815940 • 

0.184060 

1.0 

0.241971 

0.841345 

0.158655 

1.2 

0.194186 

0.884930 

0.115070 

1.4 

0.149727 

0.919243 

0.080757 

1.6 

0.110921 

0.945201 

0.054799 

1.8 

0.078950 

0.964069 

0.035931 

2.0 

0.053991 

0.977256 

0.02 2750 

2.2 

0.035475 

0.986097 

0.013903 

2* 4 

0.022395 

0.991 803 

0.008197 

2.6 

0.013583 

0.995339 

Q.OO 466 I 

2.8 

0.007915 

0.997 495 

0.002505 

3.0 

0.004432 

0,998650 

0.001350 

3.5 

0.000873 

0*999768 

0.000232 

4.0 

0.000134 

0,999968 

0.000032 

4.5 

0.000016 

0.999996 ' 

0*000004 

5.0 

0.00000.15 

0.9999997 

0.0000003 


+ Erom : Haugen, E.B, , 'Probabilistic Approach to Design* 
John Wiley and Sons, New York, 1968* 
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